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The properties of an exciton in a type II quantum dot are studied under the influence of a 
perpendicular applied magnetic field. The dot is modelled by a quantum disk with radius R, 
thickness d and the electron is confined in the disk, whereas the hole is located in the barrier. The 
exciton energy and wavefunctions are calculated using a Hartree-Fock mesh method. We distinguish 
two different regimes, namely d « 2R (the hole is located at the radial boundary of the disk) and 
d » 2R (the hole is located above and below the disk), for which angular momentum (/) transitions 
are predicted with increasing magnetic field. We also considered a system of two vertically coupled 
dots where now an extra parameter is introduced, namely the interdot distance d^. For each Ih and 
for a sufficient large magnetic field, the ground state becomes spontaneous symmetry broken in which 
the electron and the hole move towards one of the dots. This transition is induced by the Coulomb 
interaction and leads to a magnetic field induced dipole moment. No such symmetry broken ground 
states are found for a single dot (and for three vertically coupled symmetric quantum disks). For a 
system of two vertically coupled truncated cones, which is asymmetric from the start, we still find 
angular momentum transitions. For a symmetric system of three vertically coupled quantum disks, 
the system resembles for small d^ the pillar-like regime of a single dot, where the hole tends to stay 
at the radial boundary, which induces angular momentum transitions with increasing magnetic field. 
For larger d^ the hole can sit between the disks and the Ih = state remains the groundstate for 
the whole _B-region. 
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I. INTRODUCTION 

The study of self-assembled quantum dots, as realized by the Stranski-Krastanow growth mode, has been a fasci- 
nating research area during the last decade. The possibility to discover new physics in these zero-dimensional (OD) 
structures together with possible applications in opto-electronics has led to many experimental and theoretical results. 
Most studies were devoted to type I quantum dots where both electron and hole are spatially located inside 

the quantum dots. Relatively few works however were done for type II quantum dots, where electrons and holes are 
spatially separated. One expects from this type of dots interesting properties, such as large tunability of the emission 
energy and radiative lifetimes which are considerably longer than for their type I counterparts p^ . 

Within the group of type II quantum dots, one can still distinguish between two systems. Namely, the system 
which confines the electrons inside the quantum dot, and forms an antidot for the hole, like e.g. the InP/GalnP 
dots (if strain is neglected). Besides this, there are the systems, such as GaSb/GaAs and InAs/Si ||ll[], where 
the holes are confined within the dot, but where the electrons are located outside. Experimental studies on this 



type of quantum dots have been performed in Refs. [p^-|r5[ for the InP/GalnP system and in Refs. ^6 17 1 for the 
GaSb/GaAs system. Theoretical investigations on InP/GalnP type II dots were performed by Pryor et al. jl^ and 
Tadic et al. ]T9[|, who used a strain-dependent k • p Hamiltonian to calculate the electronic structure. Nomura et 
al. calculated the Landau levels in a high magnetic field, solving the Hartree equations self-consistently, while 
the present authors studied an exciton in a planar quantum dot using a Hartree-Fock approximation pT[ |. Lelong 
et al. [p2[ calculated the binding energy of excitons, charged excitons and biexcitons in GaSb/GaAs dots using the 
Hartree-Fock approximation. The magneto-exciton in a GaSb/GaAs dot was investigated by Kalameitsev et al. [ p3[ . 

In the present paper we study a single electron-hole pair bound by the Coulomb interaction, i.e. an exciton, in a 
model type II quantum dot. The quantum dot is modelled by a quantum disk of finite height, which is an extension of 
our previous work [ pl| where we did not take into account the z -direction. As an example, we consider the electron 
in the dot and the hole in the barrier material. The conclusions of our paper are also valid for the reverse situation 
where the hole is inside the dot and the electron is outside the dot. In the present work, we neglect strain effects and 
for the confinement potential we use a hard wall of finite height, confining the electron inside the dot and repelling 
the hole to the barrier region. Furthermore, a magnetic field was applied in the growth direction, i.e. B = Be^. 



1 



The first part of this paper deals with the study of the exciton in a single dot. Investigation of the influence of 
the disk parameters, namely the radius R and thickness d, showed that one can distinguish between two regimes: the 
disk-like regime and the pillar-like regime. For the first one, the disk thickness d is much smaller than its diameter, 
i.e. d << 2R, whereas the latter one describes the system with d >> 2R. Although the magnetic field couples only 
with the in-plane electron and hole coordinates, we will show that these systems, with different z-extension, behave 
very differently under an applied magnetic field. 

In the second part of the paper, we study the properties of vertically coupled quantum dots. The interest of this 
study lies in the following. It is generally known that in reality self-assembled quantum dots resemble more our 
disk-like system than our pillar-like system. However, it is also known that it is not too difficult to form vertical 
stacks of these disk-like dots and thereby create a system which could behave like a pillar-like system. 

The paper is organized as follows. In Sec. II, we give a brief discussion of our theoretical model. The numerical 
results in the absence of a magnetic field are presented in Sec. III. Section IV deals with the result for a perpendicular 
applied magnetic field. Part A of this Section discusses the results for the single dot disk-like system, whereas Part 
B is dedicated to the pillar-like system. Parts C and D deal with two vertically coupled dots, respectively disks and 
truncated cones. Parts E and F deal with three vertically coupled dots, respectively for small and large interdot 
distances. Finally, we make a small revision of the results for the single disk in Part G. Our results are summarized in 
Sec. V. In the Appendix, we discuss in more detail the numerical approach we used for the calculation of the Hartree 
integral. 



II. THEORETICAL MODEL 



We extended our previous approach which was valid for planar dots, and include the z-direction. The energies 
and wavefunctions are solved within the effective mass approximation (with nie and rrih the effective electron and 

y'^j^, <^c,e — eB/nie and LOc,h ~ eB/mh) and the Hartree-Fock (HF) single 
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where we made use of the axial symmetry by taking ^'e(re, v'ej -^e) — e*'=''''=?Ae(re, Zg) and 'i'ii{rii,(pi^,zii) = 
e^^i^'^"^ijj,^(^rh, Zh), and where the densities pe(r' , z') and ph{r',z') are given by respectively |4'e(re, </?e, ^e)!^ and 
\^h{'>'hT'^h,Zh)\^ ■ We made use of a finite difference scheme to solve the Hartree-Fock equations. More details 
about the implementation of this finite difference scheme can be found in Refs. [ p5pl| ]. Note that we consider only 
a single electron and single hole and therefore there are no exchange terms. The reason that these equations can still 
be called HF lies in the fact that the self-interaction is excluded. As confinement potentials we take hard walls of 
finite height: 
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with R the radius of the disk, and where we took Ve positive and Vh negative. Note that the only good quantum 
number is the total angular momentum in the z-direction, defined hy L = 1^ -\- Ih- 

We solved the equations self-consistently using an iterative procedure. Since only the electron is confined in the 
absence of any Coulomb interaction, we start with the free electron solution. The Hartree integrals are then integrated 
numerically as 
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where IC{x) is the complete elliptic integral of the first kind. More details about the calculation and numerical 
implementation of this integral is given in the Appendix. 

After convergence of the iteration procedure, the total energy is given by 

i^excito,, - + ^» + ^ / / "'X-'^-'^ ^-^^ (4) 

The contribution of the correlation to the total energy is neglected in HF, but for the self-assembled quantum dots, 
it is expected to be less than 2% ^ and for type II dots this will be even less. 



III. NUMERICAL RESULTS IN THE ABSENCE OF AN APPLIED MAGNETIC FIELD 

In our model system, there are two main parameters which can be varied: the radius of the disk R, and the thickness 
of the disk d. In the first part of our numerical study we will investigate the influence of these parameters on the 
exciton energy and wavefunction in the absence of a magnetic field. Following material parameters, typical for the 
InP/InGaP dot system, were taken: me = 0.077mo, nih = 0.60too, e = 12.61, = 250meV and Vh — —bOmeV. 

The variation of the disk radius R and thickness d has large consequences for the hole wavefunction. Namely, the 
hole always tends to sit as closely as possible to the electron. In this way, when the disk thickness is larger than its 
diameter {d >> 2R), the hole will prefer to sit at the radial boundary of the disk. When we consider very thin disks 
however {d << 2R), it will be much more favourable for the hole to be located above and below the disk. This state 



was not possible in our previous study [£1 
sit at the radial boundary of the dot. To i^ 
boundary of the disk as 



of the two-dimensional flat quantum dot, where the hole was forced to 
lustrate this, we calculated the probability to flnd the hole at the radial 



Pstde^2n dzh drh rh\'ifh{rh,Zh)\'^ ■ (5) 

When calculating this value for varying R and d, we obtained a phase diagram for the position of the hole wavefunction, 
which is shown in Fig. 1(a). The hole confinement potential in all these calculations was fixed to Vh = — 50meK The 
solid curve indicates where 50% of the hole is located at the radial boundary, for the dashed curves it is respectively 
25% and 75%. In order to have a more visual picture of the hole state, we present contourplots of the hole wavefunction 
\^h{'rinZh)\^ for three characteristic situations in the phase diagram, as indicated by the numbered stars. Fig. 1(b) 
shows the hole density for R = Anm, d = 12nm, where the dotted lines indicate the boundary of the quantum disk. 
We clearly see that the hole is mainly situated at the radial side of the disk. Fig. 1(c) was made for R = 8nm, 
d = 8nm and at this position, exactly 50% of the wavefunction is located at > R. The third plot. Fig. 1(d), depicts 
the result for a thin disk, with R = 12nm and d = Anm, where the hole wavefunction is mostly located above and 
below the disk, where it is nearer to the electron. 

Thus we can distinguish between two main systems: a) a disk-like system with d « 2R, where the hole will be 
located above or below the quantum disk, and b) a pillar- like system with d » 2R, where the hole sits mainly at 
the radial boundary of the disk. The solid line in Fig. 1(a) can be approximated by the curve d = 1.76i? — 5.41nm. 
We will see in the following section that a magnetic field has a very different effect on these two classes of systems. 



IV. RESULTS FOR A PERPENDICULAR APPLIED MAGNETIC FIELD 



A. Single dot: a disk-like system 



After having discussed the influence of the main parameters of our system, we are now ready to consider an extra 
feature, i.e. a magnetic fleld applied along the ^-direction, parallel to the growth direction. The magnetic field will 
squeeze the wavefunctions in the radial direction. When d << 2R, the hole is above and below the disk, and will be 
pressed stronger to the center of the disk. Fig. 2 depicts the exciton energy as a function of the magnetic field for a 
disk with R = lOnm and d = 2nm for two different values of the hole angular momentum Ih- Note that this exciton 
energy includes both the one particle energies and the Coulomb energy, i.e. Eexciton = Ee + Eh — Ecoui- In the present 
case the single hole energy Eh will be 0, as the single hole is not confined. The single electron energies at i? = OT 
and B = 50T are given by respectively Ee = 192.6TOey and E^ — 212.3meV. We find a strong enhancement of the 
exciton energy of almost 20meV for B = OT B = 50T. Due to the magnetic field, the electron and hole are both 
pushed closer to the center, which enhances the Coulomb interaction. However, also the one particle electron energy 
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increases drastically with increasing magnetic field, thus cancelling the effect of the enhanced (negative) Coulomb 
energy. The two insets show what happens with the electron (dashed curves) and hole (full curves) wavefunctions 
when a magnetic field is applied. For B = OT (Fig. 2(a)) the electron is strongly confined in the quantum disk, and 
the hole is localized above and below the disk. From the figure, it is obvious that this position is energetically most 
favourable for the hole, since the hole will be very close to the electron. Because of the very small thickness, there is 
some smearing of the electron wavefunction out of the disk in the z-direction, which attracts the hole to this position. 
The corresponding wavefunctions for B = 50T are plotted in Fig. 2(b), where again the dashed curve indicates the 
electron wavefunction, and the solid curve the hole. We see that both particles are compressed in the radial direction. 
For the electron, this leads to a further penetration of the wavefunction into the barrier. Because of this, the hole 
will be even stronger attracted to the top and bottom of the disk, which leads to a larger penetration of the hole in 
the disk. As a result, the hole wavefunction will also be slightly squeezed in the z-direction. 



B. Single dot: a pillar-like system 

The exciton energy as a function of the magnetic field is plotted in Fig. 3 for the case of i? = Anm, d = 12nm and 
Vh ~ — SOmeV. An interesting feature which appears is the occurrence of (hole) angular momentum transitions with 
increasing magnetic field. The origin of the angular momentum transitions is the following. As discussed in Sec. Ill, 
the preferred position of the hole for a disk with a thickness larger than its diameter will be at the radial boundary 
of the disk. When a magnetic field is applied along the z-direction, the hole will be pushed against the border of 
the disk. Since the disk forms a barrier for the hole, it will be energetically more favourable to jump to a higher 
angular momentum 1^ state. Notice that 1^ is an approximate quantum number. This feature was also apparent in 
our previous study of planar type II dots where we forced the hole to sit at the radial boundary. This is similar 
to an exciton in a quantum-ring structure p6| , p7| . In the present study however, we allow for the hole to relax also 
in the z-direction. Therefore, for a small thickness of the disk (disk-like system), it will be much more favourable 
for the hole to move in the z-direction and only for thick enough disks, such transitions occur. Contourplots of 
the wavefunctions are shown, respectively for the electron 9X B — OT (inset of Fig. 3) and for the hole sX B — OT, 
B = 25T, B = 50T and B = QOT (Figs. 4(a)-(d)). We find that the electron is located in the center of the disk, with 
almost no penetration into the barrier. At _B = OT, the hole is mainly sitting at the radial boundary of the disk, 
although there is some extent of the wavefunction towards the top and bottom of the disk, and there is also some 
penetration into the dot. When increasing the magnetic field, the hole jumps to the l}^ = 1 state around B ~ 15T. 
We plotted the hole wavefunction aX B = 25T (Fig. 4(b)), and we see that the hole is pushed at the side, where 
it has now more space to extend in the radial direction. Further increasing the magnetic field, leads to a jump to 

= 2 at S ~ 45T and to = 3 at B ~ TdT. The contourplots were made for respectively B = SOT (Fig. 4(c)) 
and B = 90T (Fig. 4(d)), and we see again that the hole is sitting at the side of the disk. Notice that, for all four 
plots, the probability to find the hole at r 2± 12nm is almost the same, i.e. {"^uifh = linni, z/j = 0)| ~ 0.025. In the 
absence of those Ih transitions, the magnetic field would have compressed the hole much more strongly to the radial 
boundary of the disk which leads to a larger energy as is shown clearly by the solid curve in Fig. 3. 

In order to investigate the influence of the hole confinement potential Vh on the angular momentum transitions, we 
made a phase diagram of the angular momentum state of the hole as a function of both the confinement potential 
and the magnetic field. This is shown in Fig. 5 for a disk with radius R = 6nm and thickness d = lAnm. The figure 
shows that up to Vh = —HmeV no angular momentum transitions occur. In this region Vh is too small to form a 
barrier for the hole, and the hole jumps inside the disk due to the Coulomb interaction, forming a type I system. We 
also see that the transition between type I and type II is very sharp. In our previous study of the planar dots 
we found from a similar figure a re-entrant behaviour to the Ih — state for sufficient large magnetic fields. This is 
not so evident in the present case, as we find that the line which separates the type I from the type II behaviour is 
almost vertical. 

From an experimental point of view, it is interesting to look at the probability for recombination of the exciton, 
which is proportional to the square of the overlap integral 

/flu /•QO 
*e(r)*,,(r)dr = / e''('^+''')'^d^ / i^,{r)Mr)rdr. (6) 
Jo Jo 

The integral over the angle gives 27r^i^-|_;^. This means that the probability for de-excitation is only non-zero for 
le + Ih = 0. This implies that, after an angular momentum transition, the probability for recombination decreases 
strongly, leading to a vanishing of the photoluminescence (PL) spectrum after a certain value of the magnetic field. 
Fig. 6 depicts the overlap integral / as a function of the confinement potential of the hole, for both B — OT (solid 
curve) and B = 25T (dashed curve). For B = OT the ground state will always be the Ih — state, and the overlap 
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integral will always be non-zero. However, as the system goes from type I to type II we see a strong decrease of the 
overlap, because the exciton becomes spatially indirect. For B = 25T we find a different picture: after the system has 
become type II, the overlap integral becomes zero. This is due to the fact that at i? = 25T an angular momentum 
transition has occurred and the condition + Ih = is no longer fulfilled. 

C. Two vertically coupled dots 

The inclusion of the z-direction in our calculation enables us to investigate systems of vertically coupled quantum 
dots. We have now an extra parameter to vary, namely the interdot distance dz- In this first part, we focus on two 
stacked dots of the same size, namely radius R = 6nm and thickness d = 6nm, and with dz = 3.6nm. The total stack 
height of 15.6nm is thus larger than the diameter of the disks, being 12nm. Therefore we expect to find a magnetic 
field behaviour which resembles that of the single disk pillar-like system. The result for the exciton energy as function 
of the magnetic field is shown in Fig. 7. Indeed we find angular momentum transitions with increasing magnetic field. 
The inset. Fig. 7(a), shows a contourplot of the electron wavefunction at B — OT, which is symmetrically distributed 
over the two disks. The Coulomb interaction attracts the hole as close as possible to the electron, in this case the 
radial boundary of the disks. As one can see in Fig. 7(b), the hole tends also to sit between the disks, to be even 
closer to the electron. A magnetic field acts on the wavefunctions in the radial direction, and thereby pushes the hole 
stronger between the disks. There is however not enough space, and therefore it will be energetically more favourable 
for the hole to jump to a higher Z^-state. 

When giving a closer look at the behaviour of the different Z/j-states in a magnetic field (Fig. 7), a remarkable feature 
appears: the Ih — curve exhibits a kind of kink around B — 15T, and we find a similar feature for — 1 around 
B = 45T. Investigation of the wavefunctions (Figs. 8(a-c)) learns us that we are dealing with a spontaneous symmetry 
breaking |^8|. The corresponding hole wavefunction aX B = 15T, SOT and SOT, for = is shown in Fig. 8(a,b,c) 
respectively, where one can see clearly the increase of asymmetry with magnetic field. But for such magnetic fields 
the Ih — state is not the ground state. This means that the jump to a higher angular momentum state is still 
preferred above the asymmetric state (see Fig. 8(d) for It = ^ and B = 50T). We also found that with increasing Ih 
the symmetry breaking will occur at larger magnetic fields. 

In order to investigate whether or not the ground state configuration is asymmetric, we calculated 

/•oo / /"OO A) \ 

^Pasymm = 27r / dr^ Th [ j dzh - dzh ) h{rh, Zh)\^ , (7) 

which expresses the degree of asymmetry of the hole wavefunction in the z-direction. This quantity is plotted in Fig. 9 
as function of the magnetic field for Z/i = 0, 1, 2,3. This plot confirms that after B ~ 15T the Z^ = state becomes 
highly asymmetric, and for higher B also Ih = 1 becomes asymmetric. The asymmetric states for the subsequent 1^ 
states occur for B > 50T. The contourplot of the hole wavefunction a,t B = SOT is shown in Fig. 8(d) for the ground 
state, which is Ih — 3. The inset of Fig. 9 shows the dipole moment for the different Ih states with increasing magnetic 
field. This is an interesting quantity, as it can be measured experimentally. Note that for our symmetric system, we 
have the unique feature that a magnetic field is able to induce a dipole moment! 

We can understand the spontaneous symmetry breaking as follows. In the asymmetric case the electron wavefunction 
is of course attracted to the hole (see Fig. 8(b), dashed lines) and thus stronger confined, in this way increasing its 
single particle energy. Therefore, the broken symmetry state can only occur if it can overcome this increase in energy 
due to a larger Coulomb attraction between electron and hole. We mentioned that an increasing magnetic field 
favoured the broken symmetry state. The reason is that with increasing magnetic field the electron becomes stronger 
confined in the dot, and consequently expells more in the z-direction. In the broken symmetry case it even expells 
more, resulting in a larger overlap with the hole, which can eventually make the asymmetric state lower in energy 
than the symmetric state. If this is true, than the broken symmetry would occur even sooner with increasing magnetic 
field for thinner disks, where the electron is expelled more in the barrier material in the z-direction. That this is 
indeed the case is shown in Fig. 10 where for two thin vertically coupled disks (i? = 12 nm, d = 3 nm and dz =3nm, 
Ih = 0) we find already a broken symmetry state for zero magnetic field. In this case, this broken symmetry state is 
also the ground state. 

D. Two vertically coupled truncated cones 

It is now interesting to check whether the above angular momentum transitions survive in a system which is 
completely asymmetric from the beginning, such as a system of two vertically coupled truncated cones. The two 
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cones arc equal, and defined by a base radius i?;, of Tnm, a thickness d = Gnm, and an angle of 71° (dashed lines 
in the insets of Fig. 11). The hole and electron wavefunctions (Figs. 11 (a) and (b)) show indeed the asymmetric 
behaviour, at B — OT and for — 0. The result for the exciton energy as a function of the magnetic field is depicted 
in Fig. 11. The interesting feature is that, even for this system, angular momentum transitions occur. Since there is 
still some part of the hole wavefunction located at the radial border and also some part between the disks, it appears 
to be still more favourable for the hole to jump to the higher ^;j-state than to move more and more below the disk. 
Figs. 12(a-c) show contourplots of the hole wavefunction in the ground state for subsequent i/i-states, respectively 
lh = l ai B = 20T, Ih = 2 a.t B = SOT, and = 3 at S = 50T. 

E. Three vertically coupled dots: small interdot distance dz 

For simplicity we consider only the case of three identical dots. The first system under investigation contains three 
vertically coupled quantum disks, each with disk radius R = Snm, disk thickness d = 3nm, and with interdot distance 
dz = Snm. In this case the confinement of the hole Vh was taken to be —30meV. Fig. 13 shows the result for the exciton 
energy as a function of the magnetic field. Similar to the case of the pillar-like single dot, we find angular momentum 
transitions with increasing magnetic field. The origin of these angular momentum transitions can be understood by 
looking at the wavefunctions. The insets of Fig. 13 show the electron wavefunctions at (a) B = OT and (b) B = 50T, 
and it appears that the main part is located in the middle dot, whereas there is some small extent of the wavefunction 
into both the upper and lower dots. At B = SOT, the wavefunction is more squeezed in the radial direction. The 
evolution of the hole wavefunction with increasing magnetic field is depicted in Fig. 14. At i? = OT (Fig. 14(a)), 
the main part of the hole wavefunction is situated at the radial boundary of the stacked system and its probability 
to sit at the radial side is Pside = 73%. There is some extent of the wavefunction towards the top and bottom of 
the stack, and also between the dots. When the magnetic field increases (Fig. 14(b), B = lOT), the wavefunction is 
pushed further between the disks, because it prefers to sit as closely as possible to the electron. However, there is not 
enough space between the disks to confine the whole wavefunction, and therefore it is energetically more favourable 
for the hole to jump to a higher angular momentum state. Further increasing the magnetic field leads to more angular 
momentum transitions (Figs. 14(c-f)). 

As can be seen from the wavefunctions, the ground state for this system is always symmetric. The asymmetry, as 
found in the two-disks system, is however still present, but it turns out not to be the ground state. For the Ih = 
state, the wavefunctions become asymmetric at about B = 20T, resulting in a slight bow in the full curve of Fig. 13. 
The fact that the asymmetry is less pronounced in this system can be understood as follows: in the three-disks system 
the main part of the electron is located in the middle disk. This is in contrast to the two-disks system, where the 
electron wavefunction has equal parts sitting in the two disks. Even when an asymmetry is induced in the present 
three-disks system, the electron will still be mainly located in the middle disk, and therefore the hole will sit more 
tightly at the radial boundary than was the case for the two coupled disks. 

F. Three vertically coupled dots: large interdot distance dz 

When we increase the interdot distance d^ up to 5.5nm, the hole will start to be confined between the disks. 
Fig. 15 shows that there are no angular momentum transitions for this system and it resembles very much the narrow 
disk-like system (see Fig. 2). In this figure the solid curve denotes the It = state, whereas the dashed curve denotes 
the l/i = 1 curve. We turn again to the wavefunctions to explain this behaviour. The electron is sitting mainly in 
the middle dot, this time with smaller extent into the upper and lower dots (Figs. 15(a) and (b) for respectively 
B = OT and B = 50T). The contourplot for the hole wavefunction (Fig. 16(a)) shows that already at B = OT the 
hole is completely situated between the disks. This is the preferred place for the hole, as it tends to sit as close as 
possible to the electron. Further increasing the magnetic field, squeezes the wavefunction more in the radial direction 
(Fig. 16(b)). As there is still enough space between the disks, there is no need for the hole to jump to a higher angular 
momentum state, and Ih = remains the ground state over the whole B-region. 

For this case, no asymmetry occurs up to _B = SOT. This means that this system is quite stable against spontaneous 
symmetry breaking. This stability is due to the fact that again the electron is strongly located in the middle disk, 
and that the hole is now sitting between the disks instead at the radial boundary. 
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G. Revision of the results for the single disk systems 

We now turn back to the single disk systems, both the disk-Uke and the pillar-like, where we did not talk about 
broken symmetry solutions. However also in this case they can be found. 

We find that the disk-like system, i.e. R = lOnm and d = 2nm, is perfectly stable up to S = SOT and we find 
no asymmetric behaviour. This is in perfect agreement with the result for the three coupled disks with dz = 5.5nm. 
Indeed, for the disk-like system we have again the electron sitting in the disk with the hole above and below the disks, 
which is quite similar to the second three-disks system. 

For the pillar-like system, i.e. R = Anm and d = 12nm, we do find an asymmetry, but only for the 1^ = state and 
at very high magnetic fields, i.e. starting around B ~ SOT. The difference in energy between the states with symmetric 
and antisymmetric wavefunctions is only O.OSmey for B = 90T. The ground state however remains symmetric over 
the total considered B-region. Thus the wavefunctions, as plotted in Fig. 4, remain unchanged. 

V. CONCLUSIONS 

We studied an exciton in a type II quantum disk of radius R and thickness d. For the confinement potential we took 
a hard wall of finite height, with the hole located in the barrier, which is only confined by the interaction with the 
electron, confined inside the dot. We calculated the exciton energy and wavefunction using the Hartree-Fock mesh 
method, which allows us to start without any knowledge of the single hole wavefunction. We limit ourselves to the 
study of a model system in which strain effects arc neglected. 

In the first part, we examined the case of an exciton in a single disk. First we investigated the exciton properties 
in the absence of a magnetic field. The calculation of the probability for the hole wavefunction to sit at the radial 
boundary tells us that we can distinguish two "regimes": disk-like systems with d « 2R, where the hole will prefer 
to sit above and below the quantum disk, and pillar- like systems with d » 2i?, where the hole will be located at the 
radial boundary of the disk. 

Applying a magnetic field along the 2;-direction results in a different behaviour for these two systems. For the 
disk-like system, the hole is squeezed in the radial direction by the magnetic field, and the ground state of the system 
is the Ifi = state for all values of B. For the pillar-like system, the magnetic field pushes the hole closer to the disk 
boundary, which forms a barrier for the hole. Therefore it will be energetically more favourable for the hole to jump 
to a higher angular momentum Ih state, as then the wavefunction is able to relax away from the radial boundary. 
With increasing magnetic field, we find successive 1^ transitions. 

The investigation of a system of two vertically coupled disks was done in order to compare its magnetic field 
dependence with the pillar-like single-disk system. As was expected, we find again angular momentum transitions of 
the ground state. Moreover, a new feature appeared in this study, namely a spontaneous symmetry breaking, induced 
by both the Coulomb interaction and the magnetic field. We found though that with increasing magnetic field, it is 
still preferable for the hole to jump to a higher angular momentum state. Investigation of a system of two vertically 
coupled truncated cones learned us that also in this case angular momentum transitions occur, although the system 
is highly asymmetric. 

A system of three vertically coupled disks was investigated in the third part of the paper. It was shown that the 
system with small interdot distance dz behaves similarly to that of the single pillar-like disk. Again we found angular 
momentum transitions with increasing magnetic field. For a larger interdot distance however, the hole wavefunction 
tends to sit between the disks. In this case, an increasing magnetic field does not lead to angular momentum transitions 
anymore, and the Ih = state remains the ground state. No spontaneous symmetry broken states are found as the 
ground state. This is similar to the single disk system where symmetry broken states occur only as excited states. 
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VII. APPENDIX: CALCULATION OF THE HARTREE INTEGRAL 



The Hartree integral expresses the interaction effect of one particle on another, and is given by 

^^dr' ^ fdz' f dr'r' f "^''^''^ (Al) 



y^j,2 _j_ J.I2 _ 2j-fi cos((p — ip') + (z — z')^ 



which contains a z-dependence which was absent in our previous work (see Appendix of Ref. [0). We can remove the 
(/3-dependence, since we deal with a cylindrical symmetric system. The integral over the angle becomes the complete 
elliptic integral of the first kind, which converts Eq. (Al) into 

4 / .z' / - { '-^ -\ dr'. (A2) 



The radial integral has to be solved numerically. We use a polynomial approximation for the elliptic function p9[ , 
namely 

lC{x) = [ao + aix' + 022;'^] - [&o + Wx' + 622;'^] ln(a;'), (A3) 

with x' = \ — X and where the coefficients and hi are given in Ref. p9| . The commonly used trapezoidal rule will 
give bad results, as this elliptic function implies the appearance of a logaritmic divergence in the integrand. Therefore 
we introduced the so-called 'logaritmically weighted method' which takes care of this problem. 
Generally, the following integral can be considered: 

/(r, z) = /' dxF{x) In ('f^^^l±il^^ , (A4) 
Jo \{x + rf + {z~ z'Y ) 

which, after transformation, becomes 

/(., z) ^ f dxFix + M) m f ^r'^rl^l^^r'l^ ' ^^^^ 

fr^Jo \{x + {r + hi)Y + {z - z'Y ) 

with h the discretization step and N the number of steps. If we replace F{x + hi) by Fi + — Fi] (x/h) , we can 
write (A5) as 

N-l 

I{r, z)=J2 {F^Mr, z) + [F^+l - F,]C,{r, z)}, (A6) 

1=0 

and the remaining problem is the calculation of the coefficients Ai{r,z) and Ci{r,z). The integrals which determine 
the coefficients can be solved exactly, which leads to the following results: 



with 



and 



with 



Ai{r, z) = a{h - (r - hi)) - a{h + {r + hi)), (A7) 

a(y) ^ P {y) + P {h - y) + 2\z - z'\ ( arctan , ^ „ + arctan ) , (A8(a)) 

P{y) ^y\r,(y' + {z-z'f), (A8(b)) 

Ci (r, z) = h-^[c{h -{r-hi))- c{h + (r + hi))] - 2r, (A9) 
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c(y) 



- (y i2h - yf + {z- z'f) Q (y) + - [{h - yf - {z - z'f) Q{h~y) 



+2 |z — z'l (ft, — y) I arctan -j + arctan 





(AlO(a)) 



Q{y) 



In (y2 + _ ^')2^ . 



(AlO(b)) 



Electronic mail: karenj@uia.ua.ac.be 
° Electronic mail: bpartoen@uia.ua.ac.be 
^ Electronic mail: peeters@uia.ua.ac.be 
[1] For a review, see e.g., D. Bimberg, M. Grundmann, and N.N Ledentsov, in Quantum Dot Hetero structures (John Wiley & 

Sons, Chichester, 1999). 

[2] A. Polimeni, S.T. Stoddart, M. Henini, L. Eaves, P.C. Main, K. Uchida, R.K. Hayden, and N. Miura, Physica E 2, 662 
(1998). 

[3] L.R. Wilson, D.J. Mowbray, M.S. Skolnick, M. Morifuji, M.J. Steer, LA. Larkin, and M. Hopkinson, Phys. Rev. B 57, 
R2073 (1998). 

[4] O. Stier, M. Grundmann, and D. Bimberg, Phys. Rev. B 59, 5688 (1999). 

[5] M. Brasken, M. Lindberg, D. Sundholm, and J. Olsen, Phys. Rev. B 61, 7652 (2000). 

[6] M. Bayer, A. Schmidt, A Forchel, F. Faller, T.L. Reinecke, P.A. Knipp, A. A. Dremin, and V.D. Kulakovskii, Phys. Rev. 

Lett. 74, 3439 (1995). 
[7] J. Song and S.E. UUoa, Phys. Rev. B 52, 9015 (1995). 

[8] K.L. Janssens, F.M. Peelers, and V.A. Schweigert, Phys. Rev. B 63, 205311 (2001). 
[9] M. Korkusinski and P. Hawrylak, Phys. Rev. B 63, 195311 (2001). 
[10] L. Miiller-Kirsch, R. Heitz, A. Schliwa, O. stier, and D. Bimberg, Appl. Phys. Lett. 78, 1418 (2001). 
[11] R. Heitz, N.N. Ledentsov, D. Bimberg, M.V. Maximov, A.Yu. Egorov, V.M. Ustinov, A.E. Zhukov, Zh.L Alferov, G.E. 

Cirlin, LP. Shoshnikov, N.D. Zakharov, P. Werner, and U. Gosele, Appl. Phys. Lett. 74, 1701 (1999). 
[12] C.-K. Sun, G. Wang, J.E. Bowers, B. Brar, H.-R. Blank, H. Kroemer, and M. Pilkuhn, Appl. Phys. Lett. 68, 1543 (1996). 
[13] S. Nomura, L. Samuelson, M.-E. Pistol, K. Uchida, N. Miura, T. Sugano, and Y. Aoyagi, Appl. Phys. Lett. 71, 2316 
(1997). 

[14] M. Hayne, R. Provoost, M.K. Zundel, Y.M. Manz, K. Eberl, and V.V. Moshchalkov, Phys. Rev. B 62, 10324 (2000). 
[15] M. Sugisaki, H.-W. Ren, K. Nishi, S. Sugou, T. Okuno, and Y. Masumoto, Physica B 256-258, 169 (1998). 
[16] L. Miiller-Kirsch, R. Heitz, A. Schliwa, O. Stier, D. Bimberg, H, Kirmse, and W. Neumann, Appl. Phys. Lett. 78, 1418 
(2001). 

[17] F. Hatami, M. Grundmann, N.N. Ledentsov, F. HeinrichsdorfT, R. Heitz, J. Bohrer, D. Bimberg, S.S. Ruvimov, P. Werner, 

V.M. Ustinov, P.S. Kop'ev, and Zh.L Alferov, Phys. Rev B 57, 4635 (1998). 
[18] C. Pryor, M.-E. Pistol, and L. Samuelson, Phys. Rev B 56, 10404 (1997). 
[19] M. Tadic, F.M. Peelers, and K.L. Janssens, submitted to Phys. Rev. B. 

[20] S. Nomura, L. Samuelson, C. Pryor, M.-E. Pistol, M. Stopa, K. Uchida, N. Miura, T. Sugano, and Y. Aoyagi, Phys. Rev. 
B 58, 6744 (1998). 

[21] K.L. Janssens, B. Partoens, and F.M. Peelers, Phys. Rev. B 64, 155324 (2001). 

[22] Ph. Lelong, K. Suzuki, G. Bastard, H. Sakaki, and Y. Arakawa, Physica E 7, 393 (2000). 

[23] A.B. Kalameitsev, V.M. Kovalev, and A.O. Govorov, JETP Lett. 68, 669 (1998). 

[24] M. V. Ivanov, J. Phys. B: At. Mol. Opt. Phys. 27, 4513 (1994); ibid., Phys. Rev. A 61, 22505 (2000). 

[25] F.M. Peelers and V.A. Schweigert, Phys. Rev. B 53, 1468 (1996). 

[26] S.E. Ulloa, A.O. Govorov, A.V. Kalameitsev, R. Warburton, and K. Karrai, :ond-mat/011100(: . 
[27] A.O. Govorov, A.V. Kalameitsev, R. Warburton, K. Karrai, and S.E. Ulloa, ;ond-mat/0111007 . 

[28] We allowed for broken symmetry solutions by manually pushing the electron single particle wavefunction 0.1% more to 

the positive ^-direction with respect to the negative z-direction in the first iteration step. 
[29] M. Abramowitz and L A. Stegun, Handbook of Mathematical Functions (National Bureau of Standards, Washington DG, 

1970), expression 17.3.34. 



9 



FIG. 1. (a) Phase diagram of the probability for the hole wavefunction to sit at the radial border of the disk, as a function 
of both R and d. (b), (c) and (d) are contourplots of the hole wavefunction, at the positions as indicated by the corresponding 
numbers on the main figure. The dashed lines in the contourplots indicate the boundary of the disk. 

FIG. 2. The variation of the exciton energy with increasing magnetic field, for a disk with thickness d = 2nm and radius 
R = Wnm and two values of the hole angular momentum. The insets show the electron (dashed curve) and hole (solid curve) 
wavefunction for (a) B = OT and (b) B = SOT. 

FIG. 3. Evolution of the exciton energy with increasing magnetic field, for a disk with thickness d = 12nm and radius 
R = 4nm. Increasing the magnetic field, leads to successive transitions of the hole angular momentum Ih indicated by the 
arrows. The inset depicts a countourplot of the electron wavefunction at B = OT. The dashed line indicates the disk boundary. 

FIG. 4. Contourplots of the hole wavefunction for a disk with R = 4nm and d = \2nm, at (a) B = OT, (b) B = 25T, (c) 
B = 50r and (d) B = 90T. The dashed lines indicate the disk boundary. 

FIG. 5. Phase diagram of the angular momentum transitions as a function of the hole confinement potential Vh and the 
magnetic field B, for a disk radius R = 6nm and thickness d = 14nm. 

FIG. 6. Overlap integral as a function of the confinement potential of the hole, for B = OT (solid curve) and B = 25T 
(dashed curve). 

FIG. 7. Exciton energy as function of the magnetic field, for two coupled disks with R = 6nm, d = Qnm and dz = 3.6nm. 
The insets show contourplots of the wavefunctions at S = OT, for respectively the electron (a) and the hole (b). 

FIG. 8. Contourplots of the hole wavefunction for two coupled disks, for Z/j = at respectively B = 15T(a), B = SOT (b) 
and B = 5QT (c) and for h = 3 at B = SOT (d). 

FIG. 9. Degree of asymmetry of the different Ih states as a function of the magnetic field. The inset shows the dipole moment 
as a function of the magnetic field. 

FIG. 10. Contourplots of the electron (a) and hole (b) wavefunction at B = OT for two coupled disks with R = 12nm, 
d = 3nm and dz = 3nm. 

FIG. 11. Evolution of the exciton energy with the magnetic field for two coupled truncated cones. The insets show contour- 
plots of the electron (a) and hole (b) wavefunctions at B = OT. The dashed lines indicate the disk boundary. 

FIG. 12. Contourplots of the hole wavefunction at (a) B = 2QT, (b) B = SOT and (c) B = SOT for a double dot system 
made out of truncated cones. 

FIG. 13. Exciton energy as a function of the magnetic field, for throe vertically coupled disks with intordot distance dz = 3nm. 
Increasing the magnetic field leads to angular momentum transitions, which are indicated by the arrows. The insets show 
contourplots of the electron wavefunctions at respectively (a) B = OT and (b) B = SOT. The dashed lines indicate the disk 
boundaries. 

FIG. 14. Contourplots of the hole wavefunction for (a) B = OT, (b) B = lOT, (c) B = IST, (d) B = 2ST, (e) B = SST and 
(f) B = SOT. The dashed lines indicate the disk boundaries. 
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FIG. 15. Exciton energy as a function of the magnetic field, for three vertically coupled disks with interdot distance 
dz = 5.5nm. The insets are contourplots of the electron wavefunctions for (a) B = OT and (b) B = SOT. The dashed 
lines indicate the disk boundaries. 

FIG. 16. Contourplots of the hole wavefunctions for (a) B = OT and (b) B = SOT. The dashed lines indicate the disk 
boundary. 
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